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We study the exact renormalization group (RG) in R 2 -gravity in the effective average action formalism using the 
background field method. The truncated evolution equation (where truncation is made to low-derivatives functionals 
space) for such a theory in a de Sitter background leads to a set of nonperturbative RG equations for cosmological 
and gravitational coupling constants. The gauge dependence problem is solved by working in the physical Landau- 
DeWitt gauge corresponding to gauge-fixing independent effective action. Approximate solution of nonperturbative 
RG equations reveals the appearence of antiscreening or screening behaviour of Newtonian coupling, depending on 
the higher-derivatives coupling constants. The existence of unstable UV fixed points is also mentioned. 



estartsectionsectionl0pt-3.5ex plus -lex minus - 
.2ex2.3ex plus .2exIntroduction 

It is "well-known that Einstein quantum gravity 
(QG) is not renormalizable Q. There are QG models 
which represent extensions of Einstein gravity. One 
of them, the so-called R 2 -gravity (see Q for an intro- 
duction and review) is multiplicatively renormalizable. 
However, it is most probably a non-unitary theory, at 
least in the perturbative approach. In the situation 
when consistent quantum gravity (QG) is unknown it 
is quite reasonable to study the existing gravitational 
theories as effective theories. This gives one the pos- 
sibility of estimating QG manifestations at low energy 
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scales. 

In such a way one can reduce QG to a simpler the- 
ory described by some type of scalar Lagrangian |3|, || . 
Those models are useful for describing QG in the far 
infrared domain (at large distances) || Q . 

One can consider another approach. Let us take 
non-renormalizable Einstein gravity and work with it 
as with a usual non-renormalizable effective field the- 
ory. Then the calculation of quantum corrections is 
still possible. In such a way, quantum corrections to 
the Newtonian coupling constant and to the Newto- 
nian potential have been estimated 

Finally, one can apply the exact RG |6| in the 
study of non-renormalizable theories. There has been 
recently much activity in studying different theories 
(mainly scalar ones) using the non-perturbative RG 
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(for a list of recent papers see [Q and references therein). 
Using an average effective action and the background 
field method, a non-perturbative RG study of Einstein 
quantum gravity was recently presented The RG 
equation (or evolution equation) governing the evolu- 
tion of the effective action from a scale A cut _ ff where 
theory is well-defined to smaller scales k < A cut _ ff 
has been constructed in Ref. || . Its gauge dependence 
has been investigated in |)| . It has been shown ji| that 
in the physical gauge (the Landau-DeWitt gauge) the 
Newtonian coupling constant shows antiscreening. 

In the present paper we formulate the non-pertur- 
bative RG equation (evolution equation) in higher 
derivative QG (for a review and list of references see 
[^J ) . It is widely known that such a theory, being multi- 
plicatively renormalizable and asymptotically free, has 
a perturbatively non-unitary S-matrix. Nevertheless, 
such a theory has a lot of applications. For example, 
it may lead to more or less succesful inflation |l(| . At- 
tempts to construct supersymmetric generalizations of 
R 2 -gravity have been recently made 11 , O . 



We consider higher-derivative QG as an effective 
theory, so issues of renormalizability or (non)unitarity 
are not important for us. We adopt the formalism of 
Refs. 1^, H in such a theory and construct the scale- 
dependent gravitational average action Tk [g^v] i n the 
background field formalism. A truncated evolution 
equation is obtained. We work in the physical gauge 
corresponding to a gauge-fixing independent effective 
action. Note that we make truncation of the average ef- 
fective action to the space of low derivatives functionals 
only. Even in such a simplified variant (where higher- 
derivative couplings may be considered as free param- 
eters) the calculation of nonpcrturbative RG equations 
is very complicated. 

The paper is written as follows. In the next sec- 
tion we a give very brief review of the evolution equa- 
tion (for more details see ||, ^|) and our truncation. 
Sec. 3 is devoted to the calculation of the one-loop effec- 
tive action in R 2 -gravity in the De Sitter background. 
Such an evaluation is presented in two cases: (a) a con- 
venient effective action in a one-parameter-dependent 
gauge and (b) a gauge-fixing-independent effective ac- 
tion. The results of the above calculation are used to 
obtain the effective average action in the background 
field formalism for these two cases. (In other words, we 
obtain the r.h.s. of the evolution equation in the De 
Sitter background). In Sec. 4 we perform an explicit 
truncation of the evolution equation and obtain non- 
perturbative RG equations for the gravitational and 
cosmological coupling constants. In order to avoid the 
gauge dependence problem we work there in the gauge- 
fixing independent EA formalism (see [|l3| for an intro- 
duction). The critical points of the RG equations for 
the Newtonian and cosmological couplings and running 
Newtonian coupling are discussed in Sec. 5. Finally, 
some remarks are given in the conclusion. 



estartsectionsectionl0pt-3.5ex plus -lex minus - 
.2ex2.3ex plus .2exEvolution equation for average ef- 
fective action 

We will start from a short introduction to the aver- 
age effective action approach in quantum gravity. We 
follow mainly Ref. M where more details are presented. 

The basic elements of the approach are: 

1. The background field method |^| implying that 
9nv = 9v,v + (0.1) 

where g^ v is the background metric and /i M „ is a quan- 
tum fluctuation. 

2. A scale-dependent generating functional for the 
connected Green functions 

Dh^DC^DC" exp{-S[g + h] - S gi [h;g] 
-S gh [h,C,C;g]-A k S[h,C,C;g]-S SOUIce }. (0.2) 

Let us give a description of the quantities which enter 
into Eq. (2.2). S[g + h] is the classical action of grav- 
ity theory under discussion; 5 g f denotes the gauge- 
fixing term. As we will be interested in R 2 -gravity, 
we suppose that S s i may be of the fourth order in the 
derivatives. The set of ghosts C, C includes all ghosts 
in the theory (in R 2 -gravity we have an extra ghost, 
the so-called third ghost). Finally, A k S is the infrared 
(IR) cut-off for the gravitational field and ghosts. An 
introduction to the present formalism of studying the 
average effective action has been presented in all de- 
tail in so we will not present more details of it 
here. 5 source in (2.2) is the standard action describing 
the coupling of the gravitational field and ghosts with 
the sources t^ v , , . 

Performing a Legendre transform of W k to get the 
average effective action T k [g,g], we can obtain the 
truncated evolution equation for (see || for more 
details) 

d t T k [g,g] = l -T?[(AT^[ 9l -g] +if ^l^flTK 



Y,c l TA(-MAg,g] + RfM)dtRt[g 



(0.3) 



where t = Ink, k is the nonzero momentum scale, A 
is some constant which depends on the model under 
discussion (for Einstein gravity A = k 2 ) , R k are cut- 
offs, Ci are the weights for ghosts. For the case of R 2 - 
gravity we have the Fadeev-Popov ghost with cpp = 1 
and the third ghost with weight ctg = 1/2, and, of 
course, Mpp and Mtg ar e usually known. g^ v is the 
background metric and <7 M „ = g^ v + where h^ v is 

(2) 

the quantum field. T k is the Hessian of Y k [g,g] with 
respect to g^ v at fixed g^ v . 

The next step is to specify the truncated evolution 
equation for the theory under study. We start from the 
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i? 2 -gravity (3.1) at the UV scale A cut _ ff and evolve 
it down to smaller scales k « A cut - ff ■ We use the 
truncation where the coupling constants are replaced 
by the fc-dependent functions (see (3.1)) 



— Z Nk — , 
v v 



1 ~ 1 

A -> Afc. 



(0.4) 



Note that we do not write explicitly the A: -dependence 
for the higher-derivatives coupling constants because 
we will be restricted here only to lower-derivative terms 
(i.e. up to the linear curvature term). Then, in such 
an approach (subreduction of the full set of RG equa- 
tions) , the higher-derivative coupling constants may be 
considered as free parameters of the theory. 

Then, closely following the arguments of Ref. ||, 
we get (keeping only low-derivative terms) 

T fc [g,g} = 2n 2 Z Nk J d A x^fg [-R(g) + 2\ k ] . (0.5) 

The ghost term disappears after choosing g^ ~ g^ . 
Projecting the evolution equation on the space with 
low-derivatives terms, one gets the left-hand side of 
the truncated evolution equation (2.3) as follows 

dt^k[g,g] 

= 2k 2 f d 4 x^[-R(g)d t Z Nk + 2dt(Z Nk ~X k )]. (0.6) 



The initial conditions for Z]y k ,\ k are choosen in the 
same way as in [||. The right-hand side of the trun- 
cated evolution equation (2.3) will be defined in the 
next section, following similar one-loop arguments. We 
have to note only that, unlike the Einstein gravity, the 
projectors should include the coupling constants. We 
do not give more details as they are very similar to 
those discussed in ||. 

cstartsectionsectionl0pt-3.5ex plus -lex minus - 
.2ex2.3ex plus .2exOne-loop effective action and effec- 
tive average action in R 2 -gravity 

In this section we study the one-loop effective ac- 
tion and the average effective action in higher-derivative 
quantum gravity (for a review see ^] and references 
therein). 

The classical action in Euclidean notations has the 
following form: 



S = J d A x^eR*R* + ^C^C^ 



— R 2 - 2k 2 R + 4k 2 A 
6u z 



(0.7) 



where R*R* = \e^e X/n sR^Rlp , is the 

Weyl tensor, k~ 2 = 32irG is the Newtonian coupling 
constant, e,/ 2 ,^ 2 are the gravitational coupling con- 
stants related to the higher-derivative terms in (3.1). 
It is quite well-known that the theory with the action 



(3.1) is multiplicatively renormalizable and asymptot- 
ically free (see M for a review). 

Our first purpose will be to calculate the one-loop 
effective action for the theory with action (3.1) on de 
Sitter background: 

Rp.va{i = tt: (SiiaQvp — gppgva) R, 



Rav 



12 
1 



(0.8) 



We work in the usual background field method ||, 
where the background field is given by the de Sitter 
metric, 



(0.9) 



and h^u is the quantum gravitational field. 

We will be interested in the calculation of the ef- 
fective action in a parameter-dependent gauge. Then, 
even in the case of Einstein gravity [ll| it is known 
that one has to make a change of the quantum fields: 



h = W lv 



v£ = o, 



1 



-g^vU a-, 



£u — £u 



h x a^ v 



0. 



0, 

(0.10) 



Clearly the transformation (3.4) induces a nontrivial 
Jacobian which should be taken into account in calcu- 
lating the one-loop effective action. 

The second variation of the classical action (3.1) in 
terms of the variables (3.4) is written as follows JlGj j : 

d 4 x^g-{-^h A 



3 



A 2 (m 



32j/ 2 



l -m 2 (R-AK)Uh 



(h- n<T)[ A (m 2 )A (-i?/3) 

^m 2 (R-4A)(h-na)A (0)a 
^m 2 Q (R - 4A)aA (0)A (-i?/2)a 



2 K 2 (i?-4A)e- L A 1 (-i?/4)e J 



(0.11) 



where rn 2 , = 2k 2 f 2 , m 2 , — 2k 2 v 2 . In accordance with 
O] the constrained differential operators are intro- 
duced: 

A o (X)0 - (-□ +X)<f>, 
Ai^X)^ = (-n^ + g»vX)£ v± , 
^Z (X)h^ = {-U% + 8>^X)h^. (0.12) 
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At the next step one can choose the gauge fixing term as follows: 



where 



\, = -2v,^/,;;-i(i-i0<v7 ( 



4a 2 



-□ + — 



(0.13) 



(0.14) 



and a 2 , A are gauge parameters. 

The general expression for the one-loop EA is given by 

r = S + - mdet(<5 2 ,S + S GF ) - - lndet A""' - lndet AA 1 " 

where the standard ghost operator M^ u is calculated with the help of as follows: 

1 



M„ 



, = 2i^(-D - |) + i(AT-l)V M V v | 



(0.15) 



(0.16) 



and the third ghost operator H^ v is given in (3.8). Note that when the operators , M^ v are written, the 
properties of the de Sitter background (3.2) are taken into account. 

We follow in the evaluation of the one-loop EA the results of Ref. |l(| where this calculation was performed 
in a more complicated six-parametric gauge. For our purposes, the only case a 2 — 0, A being arbitrary (a 
Landau-DeWitt type gauge), will be of interest. 

Taking into account the ghost operators and the Jacobian of the variables change (3.4), we obtain jl6| 



r« = -lndet 
2 



A ; , ( | ) A 2 ( m 2 2 + 



P + v 2 



A 



— In dct 
2 



A 



A 



A-3 



3A 2 

A (m 2 ) +m 2 



im 2 (A-4A) 



1 



- - lndet Ai — - 



R 



2 AT 2 -3 A 
1 •■- An 



R 



A 2 -3 



(4A - R) 



In det A 



R 



A- 3 



(0.17) 



Let us present Eq. (3.11) in the form of gravitational and ghost contributions: 



grav 



' ghost 



- In det 
2 



R) +lm 2 2 (R -4A) 



- In dct 
2 



An 



R 



A-3 



A (m ) 2 + m 2 



A 2 -3 

(iTTsy 



:An 



R 



A 2 - 3 



-lndet Ai 

2 



(4A - R) 



1 



lndet - lndet 



det Ai ( - ) det Ao(n) 



In 



(let A] ! ) dot A, ; 



A 



1 



R 



— + - lndet Ai — 



-lndetA (0), (0.18) 



R 



A-3 



(0.19) 



->(i) 



The gauge dependence of the one-loop effective action is clearly seen in Eqs. (3.11)- 



where r« = T^L 
(3.13). 

In order to avoid the explicit gauge dependence one can work with the gauge-fixing independent EA (for an 
introduction see [§, |l3|). An explicit calculation has been done in Ref. Jl(| with the following result: 



r (i)GFi i. lndet 

1 grav — 2 



,f 2 



i 



-lndet Ai — 



R 



Zv 2 



— In det 
2 



-A 



— m 2 (A-4A) 



1 



A 



A 



A 



A-3 
1 



A m 2 



-IndetAi — - 



(A - 4A) 



A-3 



1 



In det An (0) + - In det A 



A 



A-3 



(0.20) 



where the parameter A is fixed: A = 3/ 2 /(/ 2 + 2^ 2 ). The ghost contribution in the gauge-fixing independent 
EA formalism is given again by (3.13) with A fixed as above. Hence, we also found the one-loop gauge-fixing 
independent EA. The use of such EA solves the problem of gauge dependence of the EA (for a discussion of the 
dependence of the gauge- fixing independent EA on the metric in the space of fields in quantum gravity see [fl7|| ). 
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uantum gravity in a constant-curvature space 
131, (see [0 for a review). 



Note that the gauge-fixing independent one-loop EA in Einstein qi 
(a background like the de Sitter space) has been discussed in Ref. 

Our final goal is related to a study of the truncated evolution equation. A necessary step in such a study is 
the expansion of the average effective action in powers of the curvature. Effectively, one should use the one-loop 
effective action in such an expansion. 

However, the transition to constrained differential operators in accordance with (3.4) introduces additional zero 
modes. This leads to a wrong answer when we expand the determinants of the constrained operators in powers of 
the curvature. Therefore it is better to represent the effective action in terms of unconstrained operators. It could 
be done with the help of the following relations [O] : 



det Ay (AT) = det(-D + X) v = dct Ai(A) det A (a - ^ 
det A T (X) = det(-D + X) T = det A 2 (A) det A X (A - -^R)det A (X - 



(0.21) 



The operators from the left-hand side are unconstrained. Note also that in order to apply the relations (3.15) we 
should also rewrite the higher-derivative operators in terms of low-derivative (second order) ones. 

For simplicity, we consider below only the one-loop EA (3.14), (3.13). This EA actually describes two cases: 
for K = 1 it coincides with the standard EA (3.12), (3.13) in the gauge K = 1 and for K = 3f 2 /{f 2 + 2v 2 ) it 
describes the gauge-fixing-independent EA. 

First of all, we rewrite the ghost contribution (3.13) in terms of unconstrained operators: 



.(i) 

ghost 



l 



In det A 



v 



In 



det A 



v 



j detA ( — 



det A 



(0.22) 



Hence the ghost part is expressed in terms of unconstrained operators. 
For the gravitational part we get 



rWjf FJ =ilndet 



A, 



1 



/ b+\/b' 2 -ic 5_ 

^ 2 12 



R 



) 



In det 



- In dct 



A«, 



R 



K-3 



AoK) 



1 



fc-Vb 2 -4c 



1 . ( A v (-R/4)\ 1 , . (R 
- In det < —-7 1 — r \ + ~ In dct A v — 

2 \Ao(-i?/2)J 2 



K-3 



ml {R - 4A) 




1, , f R 
2 In dct A ( — 



(0.23) 



where 



R 



3v 2 



R 

c = — | ///- 



3l^ 2 



R 



2K 



i 2 2 (R - 4A) 



Recall that for K = 1 , (3.16) plus (3.17) gives the standard one-loop EA in the gauge K = 1. However, now this 
EA is expressed in terms of unconstrained differential operators. 

Now we can write the average effective action in the theory. First of all, to write the evolution equation we 
have to include the cut-off term Af-S . In other words, in the calculation of Wk — InZfc in the exponent of the 



path integrand we have to consider not only r^ 2 ^ av , rj, z; gh and the ghost term, but also AkS . 

The coefficients Zf rav and Z| h should be chosen so that the kinetic and cut-off terms combine to — □ + 
fc 2 ij(°)(— □ /k 2 ) for every degree of freedom. Here R^ is a dimensionless cut-off function. As in the case of pure 
Einstein gravity j^, ^| , all renormalization effects of ghosts are neglected. 

Hence the effective average action may be written in the form (see Eqs. .(3.16) and (3.17)) 



.(2) 



^k\g,9] = 2 Tr rin< z Nk 



-□ 



4c 



ilWTn \z Nk l-n + b - V ^ Tc + k^ 



\ Try In Z Nk ( -□ + b+V f^ - + k 2 R(V 



G 



L.N. Granda and S.D. Odintsov 



- Try In <^ Z 



Nk 



. D+ b-V^-c _5_ R + k2R{0) 



1 



R 



- Try In {Z Nk -□ - - + A: 2 i? (0) 



- \ Tr In |z Nfe ("□ - f + fc 2 i? (0) ) J + \ Try In jz^ ("□ + f + fc 2 i? (0) 

1 r '- D+ *±>£gE^ +fc wa 



Tr In <^ Z 



')} 



+ - Tr In < Zjvfc 



Try In 



where 
Ai = 



2 

- Tr In 
R 



-□ 



i? 



^! - y/A( - m 



k 2 R^ 



-□ + -i + k 2 R^ 
4 



i? 



Try In -□ -^+k 2 R^ 



-□ + 



i? 



K-3 



+ k 2 R^ | + Tr In | - | + fc 2 i? (0) ) } 



(0.24) 



K-3 



m, 



4mgA 



and A should be replaced by in (3.18). Thus we have got the average effective action. 

estartsectionsectionf 0pt-3.5ex plus -lex minus -.2ex2.3ex plus .2exEvolution equations for the Newtonian and 
cosmological constants 

In this section we write down the renormalization group equation (2.3) for the action (3.1). The l.h.s. of the 
truncated evolution equation is given by (2.6), where we have projected the evolution equation on the space with 
low derivatives. 

Now we want to find the r.h.s. of the evolution equation. To this end, we differentiate Eq. (3.18) with respect 
to t . Then we expand the operators in (3.18) in the curvature R because we are only interested in terms of the 
order J d 4 x^/g and / d 4 x^/gR. We also have to expand some functions of R inside the operators that appear in 
the first four terms and in the eighth and ninth terms in Eq. (3.18) up to linear terms in R. 

Let us take the first two terms in Eq. (3.18) and represent 



b±Vb 2 - Ac 



where b and c are given by (3.17). These functions may be written as 



5 + A ± y/BR 2 + CR + D 
6 



(0.25) 



(0.26) 



where A,B,C ,D depend on f ,v,A through b, c, m\ and m\. Expanding (4.2) up to terms linear in curvature 



/l,2 



R 
6 



A ± (d 1/2 + \d- 1/2 {BR 2 + CR) + ... 



A±D 1/2 + R 



\±\CD- 1 ' 2 
6 2 



(0.27) 



we can write 

fi = aiR + a 2 , 
where 



f 2 = l3iR + l3 2 



ai, 



, p 1 = 1 J 1 -±CD-V 2 )=±- 



4 3 



a 2 , p 2 = -(A±D 1 ' 2 



12 



-mi ± - 

2 2 2 



Qv 2 2 
, 8Am 



± l_(f 2 + v 2 



3v 2 



K + 6' 
6K , 





, 8Atoo " 


) m 2 


H + K \ 



-,-1/2 



K 



1/2 



(0.28) 



Here we have replaced the values of C and D . 

The functions /i and / 2 are the same for the third and fourth terms in (3.18). In a similar manner we obtain 
for the eighth and ninth terms in (3.18) the following values: 



h = J1R + 12, 



fi = 5 1 R + 6 2 
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with 



7i) <*i 



± 



2(K-3) 2(K-3) 



4 16mgA 



-1/2 



72, fa 



1 



m ± — I m 



1 / 4 16m<jA 



2A\V2 



2'"" 2 V K-S 

After this expansion we can expand the operators in (3.18) linearly in R. Let us introduce the notation 

Ar 1 (aR + A + k 2 R^) = A; 1 (x + Jfe 2 JZ<°>) - A; 2 (a + k 2 R^) aR + 0(R 2 ) (0.30) 

where a takes the values a\, (3\, 71, 5\ and A takes the values 02, 72, $2 , given by (4.4), (4.5). For a more 
compact notation let us introduce 



(0.29) 



A lX {z) = A l (X + k 2 R°{z)) 
N{z) = d t [Z Nk k 2 R^(z) ^ 

ZNk 



N (z) = d t \k 2 R^(z) 



(0.31) 
(0.32) 



Here the variable z replaces — □ /k 2 . Note that we use as a cut-off the same function as in Ref. ||: R^(z) 
z/Exp [z] — 1) . The above steps then lead to 



d t r k lg,g] = ±Tr T [JVA^J + ± Tr T 



NA?l 2 



NAyl 2 



Try [NAyl] 



NAT^l 

o~f2 



l^s [A^A^J-^Tr y [NoAyl] 



R\ iaiTr T [NA^+^Trr 



NA 



Tf3 2 



ai-^)Tr y [NAy 2 J 



Ih 



12 



NA 



V02 



±Tr s [NA; S 2 ] +i 7 iTr s 



NA" C 



+ ifc Tr s [NAs 2 2 ] + 1 Tr y [tf A^] - 



1 



1 



K-3 2 



«7: 

Tr s [AqA^q] 



(0.33) 



The terms with No are contributions of the ghosts. 

As a next step we evaluate the traces. We use the heat kernel expansion which for an arbitrary function of the 
covariant Laplacian W(D 2 ) reads 



Tt 3 [W{-D 2 )] = (47rr 2 tr,(/) |q 2 [W/] j cfixjg + ^Qi[W] J d 4 x^gR + 0(R 2 



(0.34) 



where by / we denote the unit matrix in the space of fields on which D 2 acts. Therefore tij(I) simply counts 
the number of independent degrees of freedom of the field, namely 



tr.(I) = l, £ry(/)=4, 



tr T {I) = 9. 



The sort j of fields enters into (4.6) via trj-(Z) only. Therefore we will drop the index j of A ja after evaluation 
of the traces in the heat kernel expansion. 

The functionals Q n are Mellin transforms of W: 



Qo[W] =W[0], 



Qn[W] 



r(n) Jo 



dzz n ~ x W{z) (n > 0). 



(0.35) 



Now we have to perform the heat kernel expansion (4.10) in Eq. (4.11). This leads to a polynomial in R which 
is the r.h.s. of the evolution equation. By comparison of coefficients with the l.h.s. of the evolution equation (2.6) 
we obtain in the order of J d Xy/g 



d t [Z Nk A k } 



1 1 

8^(4V)2 



5Q 2 [NA''] + 5Q 2 



■ 7Q 2 [iVAo 1 ] 



Q 2 [NA-'] + Q 2 [NAJ 2 X ] - 6Q 2 [N A~ C 



(0.36) 
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And in the order of J d A x^[g~R 
1 1 



dt(ZNk) 



24k 2 (4tt) 2 



5Qi [iVA" 1 ] + 5Qi 



P2 



7Qi [AA^ 1 ] + Qi [TV A" 1 ] 



+ Q 1 [NA£] - 12Q 1 [iVoAg 1 ] - 30 ( ai + — ) Q 2 [JVA" 2 ] - 30 ( ft + — ) Q 2 



12 



7VA7 2 

P2 



3g 2 [JVAq 2 ] - 671Q2 [iVA; 2 2 ] - 6JiQ 2 [NAjf] 



12 



if -3 



Q 2 [iV A - 2 ] 



(0.37) 



Let us introduce the cut-off-dependent integrals 



1 



* EW - r,„) J„ 



dzz 

1 

dzz 



n _ 1 R^{z)-zR^\z) 
[z + R(°){z) +w]p ' 



[z + RW(z)+w]p 

for n > 0. It follows that $o(tu) = $o(w) = (1 + w)~ p for n = 0. In addition, we use the fact that 

n . * [ggjw^ga . [2 _ t » R <o, (J) + 2oV °> 7) 



(0.38) 



(0.39) 



with ?7Ar(fc) = — dj(lnZjvfc) being the anomalous dimension of the operator ^fgR. Then we can rewrite Eqs. (4.8) 
and (4.9) in terms of <E> and $ . This leads to the following set of equations: 



d t [Z Nk A k ] 



4k 2 (4tt) 2 



fc 4 <^ 5^(a 2fc ) + 5^(/3 2fc ) - 5*2(0) + ^( 72fe ) + <$>\{5 2k ) 



dt(ZNk) — - 



- m (k) [5$J(a 2fc ) + 5*£(ft fc ) + 7^(0) + ^( 72fe ) + 



k 2 \ 10$i(a 2fc ) + 10$}(/3 2fc ) - 10*J(0) + 2$}( 72fc ) + 2$}(,5 2fe ) 

24 



(0.40) 



24k 2 (4tt) 2 

- (60ai + 5)$ 2 (a 2fc ) - (60ft + 5)$ 2 (ft fe ) + 



K-3 



6 $ 2 (0) - 12 7l $ 2 ( 72fc ) - 12<fi$2(5 2fc ) 



- i] N (k) 



5$ W) + 5$}(ft fc ) + 7^(0) + $i( 72fc ) + $J(J 2fc ) - 30(ai + ^)$ 2 (a 2fc ) 



- 30(ft + ^)^(ft fe ) - 3$ 2 (0) - 6 7l $ 2 ( 72fe ) - GS^KSik) 



(0.41) 



with 



a 2 k 


= a 2 /k 2 


— f 2K l + / 2/t fe 


(" 


P2k 


= ft/fc 2 


= Z 2 *^ — / 2,t fe 


(- 


12k 


= 72/fc 2 


2 2 , 22/ 

= V Kfe + V «fc 1 


/ 


&2k 


= <5 2 /fc 2 


2 2 2 2 / 
= 1/ K k - V K k | 


/ 



4A fc 

4A fc 
^/ 2 4 
8Afc 



1/2 



1/2 



(if-3)^ 2 
(X-3)^ 2 4 



1/2 



1/2 



(0.42) 



where k 2 , = K 2 /k 2 and A^ = A//c 2 and we have replaced the values of m 2 and m 2 , in Eqs. (4.4) and (4.5), 
respectively. Now wc introduce the dimcnsionlcss, renormalized Newtonian constant 

g k = k 2 G k = k 2 Z N \G. (0.43) 

Here G k is the renormalized Newtonian constant at scale k . The evolution equation for g k then reads 

d t g k = [2 + m (k)}g k . (0.44) 
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^From (4.18) we find the anomalous dimension r)jsr(k): 

r] N (k) = g k Bi(K k ,\ k ) + r)N(k)gkB 2 (K k , X k ) 
where 



(0.45) 



Bi(K k ,X k ) = -L(l0^(a 2fc ) + 10$109 2fc ) - 10*1(0) + 2$}( 72fe ) + 2<S>\{8 2k ) ~ (60a! + 5)$|(a 2 k) 



(60/3i + 5)*2(#j fc ) 



24 



if -3 



- 6 $^(0) - 12 7l $^( 72fc ) - 12^$^(5 2fe ) 



S 2 ( Kfe , A fe ) - -^{s^M + 5$}(/3 2fe ) + 7$J(0) + $i( 72fe ) + ^(M - 30(a x + ^)^(a 2fe ) 
- 30(/3i + &)*£(&*) - 3*1(0) - 6 7 i$l( 72fe ) - Q^ 2 {8 2k ) 



(0.46) 



g k Bi{n k ,X k ) 



Solving (4.21), 

1 - g k B 2 {n k , Afe) 

we see that the anomalous dimension t^jv is a nonpertubative quantity. From (4.16) we obtain the evolution 
equation for the cosmological constant 



(0.47) 



8 t {X k ) = - [2 - m (k)} Afe + g- jl0$ 2 (a 2fe ) + 10^G3 2fc ) - 10$ 2 (0) + 2$ 2 ( 72 fe) + 2$ 2 (<5 : 
- m (k) \b^\{a 2k ) + 5$ 2 (/3 2 fe) + 7^(0) + $ 2 ( 72 fe) + $l(5 2k ) 



2k I 



(0.48) 

Eqs. (4.20) and (4.24) with (4.23) determine the value of the running Newtonian constant and cosmological constant 
at the scale k << A cut _ ff. The above evolution equations include nonperturbative effects which go beyond a 
simple one-loop calculation. This is particularly obvious if one expands (4.23) for small values of g k : 



r) N = g k Bi{n k ,X k ) 1 + g k B 2 (K k ,X k ) + glB 2 (n k ,X k ) + 



(0.49) 



estartsectionsectionl0pt-3.5ex plus -lex minus -.2ex2.3ex plus .2exCritical points and the running Newtonian 
coupling constant 

In the present section we give some remarks about the properties of nonperturbative RG equations. First of 
all, let us estimate the qualitative behaviour of the running gravitational coupling constant. 
The dimensional quantity G k evolves according to 



d t G k = rj N G k 



(0.50) 



The set of RG equations for the coupling constants is too complicated and cannot be solved analytically. Hence 
we assume that the cosmological constant is much smaller than the IR cut-off scale, Afe << k 2 , so we can put 
Afe ~ 0. This simplifies Eqs. (4.29), (4.24) and we are left with only Eq. (5.1). After that we perform an expansion 
in powers of (G 2 )~ 1 keeping only the first term (i.e. we evaluate the functions $£(0) and $^(0)) and finally 
obtain (with g k ~ k 2 G) 

G fe = G [l-wGk 2 + ...] 

where 



(0.51) 



» = 4*1(0,0) = ^ 



50 + 22 



7t^_ 
3 



In obtaining w we use the same cut-off function as in Ref. ||. 

For Einstein gravity in the same formalism (also using the gauge-fixing independent EA) we have got [9] : 



36 



108 



1 



(0.52) 

In the case under discussion we see that the sign of w depends on the higher-derivative coupling constants: 

7tt 2 22 f 2 

w>0 if 50- — + ^->0. (0.53) 
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The coupling constant v 2 maybe chosen to be nega- 
tive (see M). So, e.g., for f 2 — 1, v 2 — ±1 we get 
w > and the Newtonian coupling decreases as k 2 
increases. In other words, we find an antiscreening be- 
haviour of the gravitational coupling. On the contrary, 
for f 2 = 1, i/ 2 = -1/2 we get w < and a screening 
behaviour for the Newtonian coupling (for the one-loop 
behaviour of the Newtonian coupling in R 2 -gravity 
with matter, see [^9|). This means that in such a 
phase the gravitational charge (mass) is screened by 
quantum fluctuations, or, in other words, the Newto- 
nian coupling is smaller at smaller distances. The sign 
of a quantum correction to the Newtonian potential 
will be different too. Note that our solution (5.2) is 
actually qualitative, and the full RG system should be 
analyzed for a better result. 

Our main qualitative result is that R 2 -gravity con- 
sidered as an effective theory may change the low- 
energy gravitational phenomena as compared with Ein- 
stein gravity. 

Let us now investigate the problem of existence of 
critical points in the theory under study. We search 
for points at which the r.h.s. of Eqs. (4.20) and (4.24) 
are equal to zero (supposing = k 2 G). For Einstein 
gravity such a study may be carried out quite easily 
(since the functions B\ and B2 depend only on \k ) ■ 
A numerical analysis of the corresponding RG system 
which is written in the physical Landau-DeWitt gauge 
gives Ref. 0: 



A fc = 0.352, 



g k = 0.348. 



(0.54) 



These points actually correspond to UV-stable fixed 
points. Note that the solutions (5.5) do not give a so- 
lution to the cosmological constant problem as a result 
of the non-perturbative RG behaviour. 

In R 2 -gravity the situation is much more compli- 
cated because the functions £>i and B 2 depend on 
Kfe and Xk and because there are the higher-derivative 
coupling constants as free parameters of the theory. 
Supposing gk — k 2 G, we can get the following unsta- 
ble fixed points for f 2 = 1/10 and v 2 = —1/4: 



A fe = 4.47, 



9k = 4.46. 



(0.55) 



Note that for other values of the higher-derivative 
coupling constants one can get numerically other val- 
ues for the unstable fixed points. 

estartsectionsectionl0pt-3.5ex plus -lex minus - 
.2ex2.3ex plus .2exDiscussion 

In the present work we have studied the truncated 
evolution equation in higher-derivative quantum grav- 
ity. Making a truncation to the space of low-derivative 
functionals, we have obtained nonperturbative RG 
equations for the Newtonian and cosmological cou- 
pling constants. The necessary step in such a study is 
the calculation of the effective average action on some 
background (we have used the de Sitter background). 



The properties of the nonperturbative RG equa- 
tions (like the existence of critical points and the be- 
haviour of the Newtonian coupling) are discussed. As 
we showed, the higher-derivative QG may behave at 
low energies qualitatively different from Einstein QG. 

The next open problem in such an approach is to 
make a better truncation of the evolution equation to 
the space of functionals with higher derivatives. In 
such a way one could obtain a complete set of nonper- 
turbative RG equations for all coupling constants: f 2 , 
v 2 , k 2 , A. Hence, unlike the present study where f 2 
and v 2 are free parameters, we might define the critical 
points of the complete phase space (the RG equations 



for 



and A are, of course, the same) 



However, in order to find the effective average ac- 
tion with a truncation to the space of higher-derivative 
functionals, we have to perform a calculation of the 
one-loop effective action in a background where R 2 
and C 2 a p may be distinguished. Clearly the de Sit- 
ter space does not belong to this class of backgrounds. 

As far as we know (see for a review) , the one-loop 
effective action for R 2 -gravity has been found only in 
the de Sitter or flat backgrounds. Even such a calcula- 
tion is extremely complicated. A generalization of such 
a result to a more complicated background (say, of the 
above sort) being, in principle, possible, is extremely 
complicated. Moreover, that is just one step in writing 
the r.h.s. of the evolution equation. After that, much 
more work is required to obtain explicitly the nonper- 
turbative RG equations for v 2 and f 2 . Hence, this 
problem is left for future research. 

Another related problem is the gauge dependence 
of the average effective action. In order to solve this 
problem, one has to do an even better truncation which 
includes all gauge parameters as independent functions 
of k. Hence, in addition to the four RG equations 
for the coupling constants, one should write some RG 
equations for all gauge parameters. Then the RG equa- 
tions for the gauge parameters should lead to some sta- 
ble fixed points. These fixed point values for the gauge 
parameters should be used in the RG equations for the 
coupling constants. It is clear that such a programme 
is too complicated and cannot be realized. 

However, there is a simpler way which we have ac- 
tually used in this work. In a study of the effective 
average action for the Yang- Mills theory (see |20 ) 
a k -dependent gauge parameter was used. It has 
been shown that there exists an attracting fixed point 
of the truncated evolution for the gauge parameter. 
This fixed point corresponds to the so-called Landau- 
DeWitt gauge. Hence, the alternative easy way of 
studying the truncated evolution equation is to work 
in the physical Landau-DeWitt gauge (actually it cor- 
responds to a study in the gauge-fixing-independent 
effective action formalism) . Similarly in Einstein grav- 
ity, in order to avoid the introduction of a tedious 
additional RG equation for the gauge parameter, one 
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can work in the Landau-DeWitt gauge which again 
corresponds to the gauge-fixing-independent EA for- 
malism (see p). In the same way, we have used here 
the gauge-fixing-independent EA in order to solve the 
gauge dependence problem for R 2 -gravity in our for- 
mulation. 

Hence, our study which indicates the qualitative 
difference between Einstein and R 2 -gravity even at 
low energy scales is a necessary step in the formula- 
tion of better truncations of the evolution equation in 
R 2 -gravity. Moreover, it is expected to be useful also 
in the studies of supersymmetric R 2 -gravity in a non- 
perturbative approach. 

estartsectionsubsection20pt-3ex plus -lex minus - 
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